Engineering directed excitonic energy transfer 
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We provide an intuitive platform for engineering exciton transfer dynamics. We show that careful 
consideration of the spectral density, which describes the system-bath interaction, leads to opportu- 
nities to engineer the transfer of an exciton. Since excitons in nanostructures are proposed for use in 
quantum information processing and artificial photosynthetic designs, our approach paves the way 
for engineering a wide range of desired exciton dynamics. We carefully describe the validity of the 
model and use experimentally relevant material parameters to show counter-intuitive examples of a 
directed exciton transfer in a linear chain of quantum dots. 
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The widely-applied Forster theory for energy transfer 
links experimental results to estimates of system infor- 
mation, particularly in biological and nanoscale appli- 
cations [ll, [^j. The usefulness of this theory is partly 
due to the simple expression of the kinetic rate constants 
as a product of electronic coupling and a spectral over- 
lap factor which captures the complexity of the environ- 
ment. Forster theory describes transport in the inco- 
herent limit, but a complementary and more elaborate 
approach, such as Redfield theory, is often required to 
describe energy transfer. However, the information es- 
sential to understanding the dynamics is buried within 
the structure of the equations. In this letter, we employ 
a quantum kinetic rate approach to distill the informa- 
tion contained in equations into a simple, yet instructive, 
formula. We use this approach to design directed exciton 
transfer mediated by an environment. 

Excitonic energy transfer (EET) has been studied in 
systems as varied as quantum dot (QD) nanostructures 
polymer chains [5|, and photosynthetic complexes 
Many applications of EET would benefit from 
controlling exciton dynamics. Perfect state transfer, as 
studied in the quantum computing community, is achiev- 
able in certain engineered systems, but only at particular 
times during coherent evolution Recent works have 
shown that environment-induced decoherence can alter 
exciton dynamics [9l-[ll|. although controlling the trans- 
fer direction has only been achieved using external poten- 
tials [l2|. Our paper builds upon the idea of engineering 
exciton transfe r by designing appropriate system-bath in- 
teractions [H, [iJI- We show that it is possible to design 
experimentally realizable systems where the environment 
can be used to direct the flow of energy. 

The Hamiltonian used in our simulation aims to cap- 
ture dynamics in a single-cxciton manifold fl5| interact- 
ing with an environment. 
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This representation is in the site basis {|s„)} of local- 
ized excitations on each of N sites, (e.g. QDs or chro- 



mophores), with excitation energy En for each site and 
inter-site coupling Jmn- The environment is described by 
a phonon bath, 



7?fa = ^fic^q(6t6q-Hl/2), 
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where 5^ (&q) is the creation (destruction) operator for a 
phonon with wavevector q. The system-bath interaction 
is assumed to be linear. 
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where describes the site-specific coupling of electronic 
and vibrational degrees of freedom. We ignore the off- 
diagonal terms in the above equation, which correspond 
to phonon- induced modulations of the inter-site coupling, 
Jmn [Hi ■ While fluctuations of the gap give rise to the di- 
agonal electron-phonon coupling considered here [l6j , the 
inter-site couplings are usually one or two orders of mag- 
nitude smaller than the excitation gap and are therefore 
kept constant. In general, the validity of this approxima- 
tion is still an open question and its applicability varies 
from system to system. 

To describe the excitonic quantum dynamics we use 
Redfield theory p7l - [l9| . which is a reduced density ma- 
trix approach in the regime of weak system-bath cou- 
pling. The formalism involves second-order perturbation 
theory in the system-bath interaction, Hgb- This method 
assumes the Markov approximation, no initial correla- 
tions between system and bath degrees of freedom, and 
a thcrmalizcd bath. We avoid the frequently employed 
secular approximation (Bloch equations) [3, [l^; it is 
important to note that for systems where the time scale 
\u!ab — ^cd\~^ is comparable or larger than the charac- 
teristic decoherence time, the coherence to population 
transfers contribute considerably to the dynamics of the 
system and must therefore be included. 

The equation of motion for the density operator p(t) 
in the excitonic energ y b asis representation, \ea) ~ 
ea \ea), is given by [Ulil, 
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with pab{t) = {ea\p{t) I^f)) and Uab = (eq - eb)/^- The 
first term on the right hand side of Eq. [5] describes the 
fully coherent dynamics in the absence of Hsb and the 
second term describes the irreversible dynamics from the 
interaction with the phonon bath. 

Correlations in bath-fluctuations on different sites [l^l 
are taken into account by using the relation g^g^* = 
3q 5q* = g^g-^^^/flcorr^ where Rmn = |Rm«| is the dis- 
tance between the sites and -Rcorr is the phonon correla- 
tion length HH . 

Using this relation for the electron-phonon couplings, 
the cross-correlation C„m(w) can now be written as, 

CrrM = e-^-/^— (6) 

where the frequency correlation function C{uj) — 
27r[n(w) + 1]( J(w) — J(— cj)), with n(w) the Bose-Einstein 
distribution and the spectral density of the bath J{uj) = 
^q|5q|^Wq(5(w — Wq), whcrc S{uj) IS thc Dirac delta func- 
tion. 

For the Hamiltonian specified in Eq. [2 the Redfield 
tensor elements are given by 

Rab,cd = ^db,ac{^^ca) + ^caMi^db) 

- Sbd ^ TaeM'^ce) - Sac ^ T^be.edi'^de), '''^^ 

e e 

where Sij is the Kronecker delta and 

^ab,cd{i^dc) = ■^Cab,cdC {cddc) 

27r \J-oo ^dc-i^ J 

Cab,cd = J2^U-')anUnbiU-^)cmUmde-'^"-^'^-', (9) 
n,m 

with V denoting the Cauchy principal value of the in- 
tegral and Una = (snlca) the transformation matrix el- 
ements relating the site basis {|s„)} and the excitonic 
basis {|ea)}. Since \ea) = J2n^ria |s„), then \Una\'"^ can 
be interpreted as thc contribution of the n-th site to the 
a-th eigenstate of Hg- 

Hereafter, we will focus on QDs at low temperature 
(10 K) as our prototypical experimental realization. We 
neglect possible inversion asymmetry of the crystal, and 
therefore the contribution of the piezoelectric coupling, 
and focus instead on the deformation potential coupling. 
As shown by Calarco et al. [l^j, the spectral density 
describing this coupling (in thc absence of an external 
electric field) is given by 

J(tj) = e(w)?7w3e-"'/"', (10) 

where, 8(a;) is the Heaviside step function. We use 
typical values for GaAs QDs for our numerical simula- 
tions [H [Hi, giving ?7 = ''i^ipu'^n " ^-^^^ 



ujc = y/2vFJP' = 1.41 ps""'^, where Df, (Dh) is the de- 
formation coupling potential for electron (hole), u the 
speed of sound within the quantum dot, p its mass den- 
sity and I the ground state localization length, assumed 
to be the same for electron and holes. We also assume a 
correlation length, i?corr, of 3 nm. 

Using Eqns.[7][ni thc population transfer rates, kab-, be- 
tween the eigcnstates a — > 6 is given by, 

kab ~ Rbb.aa ~ CabMa^i^ab)- (H) 

This equation is central to our insight into designing exci- 
tonic transfer. Though derived for a different regime, the 
form of Eq.[Tl]is similar to the widely used rate equation 
in Forstcr theory (incoherent limit), fc^^o»'«*er ^ J^^^Iab, 
where J is the electronic coupling, and / is the spectral 
overlap integral [25j . Both equations are a product of two 
terms: one dealing primarily with the description of the 
system (C in this paper or for Forstcr theory), and an- 
other largely depending on and arising from thc system- 
bath interaction (C(aj) in this paper or / for Forstcr the- 
ory) . Forster theory has been applied in fluorescence res- 
onance energy transfer (FRET) to design chromophores 
for biosensing assays 0, [2^, [131 and was recently verified 
experimentally for semiconductor QDs [1^. We use the 
simple structure of Eq.[TT]to gain microscopic and exper- 
imentally relevant insight into engineering directed and 
optimized EET. 

While C{LLiab) depends on the overlap of system 
eigenenergies with the spectral properties of the phonon 
bath (e.g., lattice vibrations, solvent, protein environ- 
ment, etc.), C depends on the transformation matrix U in 
Eq. [9l determined by the relative magnitude of electronic 
couplings with respect to site energies, site connectivity, 
and spatial correlation between sites. The aim is then to 
maximize (minimize) the product of these two factors in 
Eq. [TT]to favor (suppress) the desired rates. 

To illustrate thc applicability of Eq. [TTl we choose two 
three-site examples to highlight the importance of the 
phonon bath interaction to achieve directed EET. Since 
multiplying a system Hamiltonian by a scalar does not 
change the maximum exciton transfer probability in the 
coherent limit, thc distances and site energies for the two 
cases are chosen so that the Hamiltonians arc related 
by a multiplicative factor (3.5/f^ = H"^)- As a conse- 
quence, the C values are roughly the same for cases 1 
and 2 (Table slight differences are introduced by the 
bath correlation term in Eq. [9l In the cases chosen, the 
fully coherent evolution gives a maximum probability of 
finding an excitation of 5% (1%) for site 1(2). Therefore, 
any difference between the two examples is a result of 
interaction with the environment. 

In contrast to thc ^ factors, the C(wat) values for cases 
1 and 2 at 10 K differ by at least one order of magnitude 
due the position of the transition frequencies with respect 
to the spectral density (Table [J). Any changes to the 
system Hamiltonian can affect C(wab)- In our examples, 
the scalar factor mentioned above changes the C{ijJab) 
such that the largest population transfer rate is switched 
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from site 1 to site 2 since Wg^ ~ wfj and C{ujIi) ~ C{uj^2) 

(Fig.HD. 

TABLE I: Contributions of the system factor, C^ab,ba, and of 
the overlap between transition frequency and phonon bath 
spectral properties, C{ijJab), to the calculation of the quantum 
kinetic rates kab from energy eigenstate \ea) to \eb). The two 
cases considered are described in Fig. [T] 







Case 1 


Case 2 




-> let) 


log ( log C log fc 


log log C log k 


3 - 


^ 1 


-1.5 11.6 10.0 


-1.5 7.9 6.3 


3 - 


-> 2 


-2.4 10.7 8.3 


-2.4 11.6 9.1 


2 - 


■> 1 


-4.3 11.4 7.1 


-4.3 10.2 5.8 



From the values of ^ and C{u!ab), it is clear that case 1 
is designed such that an excitation starting on site 3 will 
tend to transfer to site 1, but in case 2 the population will 
go to site 2, albeit at different rates. Simulations of the 
quantum dynamics according to Redfield theory confirms 
this result (Fig.[2|). Over the typical exciton lifetime of 1 
ns in QDs, we not only achieve directed transfer, but also 
population enhancement compared to both the maximum 
site population during fully coherent dynamics and the 
population expected at thermal equilibrium. 

Moreover, while there are always experimental limi- 
tations in tuning parameters, the structure of Eq. [TT] is 
valuable since it partitions the effects due to the system 
and bath. Using the calculated rates and/or a visual in- 
spection of C{u!) (Fig. IB), it is easy to determine the 
impact of varying a system parameter on exciton trans- 
fer. Future work will address exciton and electron trans- 
fer between sites with varying spectral density functions, 
as well as the role of aligning dipole moment orienta- 
tions in engineering EET. We are also working to identify 
regimes in which preserved coherences enhance or reduce 
the efficiency of excitonic transfer. Of course, in situa- 
tions where multiple excitons are present in the system 
due to incident light intensity, frequency range, and/or 
optical spectral density of the quantum dots, the Hamil- 
tonian used to describe the system needs to be expanded 



accordingly (29[. To the extent that each exciton cou- 
ples to the environment through a spectral density as 
described in this letter, some of the intuition developed 
here should be transferable to these systems. However, 
there are a number of interesting complications, including 
many-body interactions among the excitons and correc- 
tions to the rates due to multi-phonon processes. A more 
careful analysis is needed when the offset in QD excita- 
tion energies is large compared with the characteristic 
frequency expanded by the phonon spectral densities or 
for systems at much higher temperature, e.g., room tem- 
perature, where multi-phonon processes are expected to 
be relevant, or even dominant. These effects are cur- 
rently under investigation and are beyond the scope of 
the current communication. 

In summary, we develop a framework for engineering 
environment-assisted and directed excitonic transfer in 
a network of coupled QDs based on a quantum kinetic 
rate approach. We emphasize the importance of how 
characteristic frequencies of a system fit within the spec- 
tral bath structure. Our examples utilize the factored 
and intuitive form of the population transfer rates equa- 
tion, which separates the contributions from the system 
(electronic) and bath (vibrational) degrees of freedom. 
This equation is similar in spirit to the rate equation for 
FRET, making it convenient to design interesting scenar- 
ios for environment-assisted transfer. Although we focus 
on QD examples, the principles presented here form the 
basis for engineering a wide range of desired EET in a 
variety of nanostructures or artificial molecular photo- 
synthetic units. 
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FIG. 1: Two cases of a three-site system are considered, a) Scaled schematic of system spacing and energy levels with details of variable site energies (En) and intersite 
Forster coupling strengths (Jmn)- For QDs with transition dipole moments aligned perpendicular to Rmn, Jmn = 100meV/-R^„, with Rmn in nanometers [sO]. b) 
Frequency correlation function for a super-ohmic spectral density. Energy basis transition frequencies for case 1 (2) are indicated by solid (dashed) vertical lines. 




FIG. 2: Site basis population probabilities for an excition starting on site 3 for case 1 (a) and case 2 (b) demonstrate the change in transfer dynamics obtained by 
scaling the Hamiltonian (simulation at 10 K). Dashed lines indicate site populations at thermal equilibrium. Energy basis coherences for case 1 (c) and case 2 (d) 
display characteristic oscillations and damping over the course of 1 ns, a typical recombination time in QD systems. 
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